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£SJ . A study of inclusive production of single hadrons with finite transverse momentum 

in neutral current deep-inelastic scattering has been carried out. Cross sections have 
been calculated using perturbative Quantum Chromodynamics at next-to-leading order 
and compared to HERA data. Predictions for light charged hadron production data 
were also calculated for large values of Q 2 to show the possible effects of the Z-boson 



i exchange. 

^ <—| 1 Introduction 

JL , The perturbative QCD approach for computing hadronic cross sections is based on the par- 

ton model picture, in which the cross section for any hard scattering process can be written 
\ as a convolution of structure (f a (x, Q 2 )) and fragmentation (D a (x, Q 2 )) functions of partons 

(quarks and gluons) and a hard-cross section factor [T]. The structure and fragmentation 
functions are non-perturbativc, universal quantities, while the hard-cross section can be cal- 
■ culated within perturbative QCD to the lowest order in the running coupling a s {Q) as long 

as Q ^ A, where A is the QCD scale. 

This naive parton model corresponds to the leading order (LO) approximation. However, 
due to the perturbative nature of a s , the running of the coupling constant could be hidden 
in higher order corrections and therefore, the LO calculation can only predict the order of 
magnitude of a given cross section. The accuracy of the perturbative QCD expansion is then 
controlled by the size of the higher order contributions. Any perturbative QCD prediction 
| needs then, next-to-leading (NLO) corrections and NLO definitions of the running coupling 

constant, and the structure and fragmentation functions. 

In this work we have carried out a NLO calculation of single hadron production in deep- 
inelastic scattering (DIS), using standard PDF sets and fragmentation functions fitted from 



H ' e + e data. 
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2 Theoretical formalism 

We are interested in single production of hadrons in DIS: 

l(k)+p(P)^l(k') + h(p h ) + X, (1) 

where the lepton can be an electron, muon or an (anti)neutrino, and the exchanged vector 
boson a photon, or Z. In this work we concentrate only on the neutral current case. 
The cross section can be written as a convolution of the partonic cross section with the 
parton distribution functions and the fragmentation functions [2] : 

da h f 1 dx fUz y f p( x 2\ da ab D h (* »\ (2) 



dxdydzdcf) J x x J z z \x ) dxdydzd(f> \z 
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where we have used the kinematic variables 

X 2P-q' V P.Jfe' Z ~ P-q 

Q 2 - Pa ■ Pb 

x = 7, . y = y, z= . 3) 

2p a • q Pa' q 

The partonic cross section may be split into leptonic and hadronic parts, 

da ab a 2 v 

V -X ab L^Hf v , (4) 



dxdydz 16n 2 Q 



where 



2y V y / s r t V 2/ / s ft V y J 

(5) 

with Sh — The coefficients A a (, contain the information about how the fermions couple 

xy 

to the vector bosons: 

X T ab L = e}-2e f v f v eXz (Q 2 ) + (a 2 f + v})(al + v 2 e )x 2 z (Q 2 ) (6) 
\£ b = -2e f a f a e xz(Q 2 ) + ia f a e v f v eX 2 z(Q 2 ), (7) 

where 

Oe,/ = T e 3 /; (8) 
« eJ = T e 3 / — 2e e j sin 2 6W, (9) 
1 Q 2 

* Z (Q 2 ) = a ■ 2 a YT~ r>2 1 ■ ( 10 ) 

4 sin 6% cos t'vK Q + M z 

T 3 y is the weak isospin and 9w the weak mixing angle. 
3 Comparisons with HERA data 

In this section we present our numerical results for single hadron inclusive production mea- 
surements from HERA data. Cross sections are calculated to NLO in the MS scheme. We 
set the number of active quark flavours nj = 5. We used the CTEQ6.6M PDF set of rcf. [3j 
and their value Aq^ d = 226 MeV. 

At first order in as, the only partonic subprocess contributing to LO is 7*, Z a + q — > q. In 
order to calculate the cross section at NLO we must include virtual and real corrections to 
this process, and using dimensional regularization and standard renormalization techniques, 
cancel all singularities that there appear. In this case the factorization/renormalization scale 
was chosen as fi = Q. 

New charged particle production data was released in 2007 by the HI collaboration |3j, that 
we can use for comparison. There, the cross section differential in the scaled momentum 
x p = -^k, normalized to the total cross section, was measured for different bins in x p and 
for 100 < Q 2 < 10000. Our results arc shown in Figure [U where the solid central line is our 
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main prediction, including the Z boson contribution, and the dashed line is the prediction 
with only virtual photon exchange. We have used the AKK set of fragmentation functions 
to make our predictions [5] . The two solid lines at the top and bottom of the central line are 
the predictions obtained when the scale was changed to fi = \\x and /z = 2/z respectively, 
which gives us an idea of the theoretical error in our calculation. 



Unfortunately, although for these val- 
ues of Q 2 , wc should be able to 
see an effect due to the presence of 
the Z boson, since these cross sec- 
tions are normalized to the total cross 
section, and the effect is present in 
both, they cancel each other out and 
we can not distinguish between the 
curve with and without the Z bo- 
son contribution. What we can do 
is plot the ratio of the unnormalizcd 
cross sections and there we should be 
able to sec what kind of effect the 
Z boson introduces. The ratio for 
the total and the differential cross sec- 
tions is shown in Figure [2 The ef- 
fect of the Z boson contribution is 
found to be of about 15% in this 
case. 
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Figure 1: Comparison of theoretical predic- 
tions using AKK with HI data. 



At second order in a S) the LO contributions 
are 



l\Z a + q- 
l\Z a + q- 



q + g, 
g + q, 
> q + q. 



(ii) 
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and we should again calculate virtual and 
real corrections to these processes. In this 
case, the cancellation is not straightfor- 
ward as in the previous case and a suit- 
able method must be used in order to cancel 
all singularities. We have chosen the dipolc 
subtraction method pQ. After carrying out 
the renormalization procedure for the vir- 
tual corrections, all the ultraviolet singular- 
ities should cancel and we have only soft 
and collinear divergences left. The idea of 
the subtraction method is to use the identity 



Figure 2: Ratio between cross sections with 
and without Z boson contribution to first or- 
der in a.,. 



NLO 



{da 1 



da J 



m+1 



da J 



m+1 



da" 



(12) 
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a NLO 



(13) 



where the approximated cross section da A is such that is has the same singular behaviour 
in d dimensions as do R . In this way, da A acts effectively as a counter term for da R and 
introducing the phase space integration, we can safely perform the limit e — > under the 
integral sign in the first term on the right-hand side of equation[T2land therefore, perform the 
integral in four dimensions. All the remaining singularities are at this point associated with 
the last two terms in equation 1 121 If we then integrate da A analytically over the one parton 
subspacc leading to the e poles, we can combine these poles with the ones coming from the 
virtual corrections and in this way, cancel all the singularities. Then we can perform the 
limit e — > and integrate numerically the rest over the m-parton phase space. The final 
structure of the calculation is 

/ [Me=0-(^)e=o]+/ W+f*** 
J m+1 J m L ^ 1 

Using this formalism, we can then obtain predictions for single hadron production at NLO at 
second order in a s . However, since we are interested in perturbative QCD effects, we require 
the hadron to carry non-zero transverse momentum (p^) in the centre-of-mass frame of the 
vector boson and the incoming proton. Two measurements are useful for comparison, both 
from the HI collaboration ( [6] and [7]). The first one is tt° production, where we have used 
again the AKK set of fragmentation functions for our predictions and p^ > 2.5 GeV. For the 
second one, D* ± production, we have used the KKKS set of fragmentation functions [8] and 
p* T > 2.0 GeV. More details about the kinematic region used in each of these measurements 
can be found in the experimental papers respectively. Our results are shown in Figures [3] 
and [4] Here again the solid red line represents our prediction and the dashed lines at the top 
and bottom of it are the corresponding predictions with a different scale as in the previous 
case. It is not possible to appreciate any effect of the Z boson contribution, since Q 2 only 
goes up to 70 GeV 2 . However, we can imagine we extend the range of Q 2 for these two 
cases, and as before, we calculate the ratio between the cross sections with and without the 
Z boson. Our results are shown in Figure [5l where we can see contributions of up to 7% for 
Q 2 > 10000. 



4 Conclusions 



The NLO cross section for single hadron production was calculated in perturbative QCD for 
neutral current DIS using the dipolc subtraction method. Results agree with the data for ir° 
and D*^ production using AKK and KKKS fragmentation functions. Effect of Z° boson in 
the cross sections found to be up to 15% in the first order calculation and 7% second order 
calculation. 
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Figure 3: Cross section differential in Q 2 and xb for 7r° production. 
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Figure 4: Cross section differential in Q and xb for D** production 



Figure 5: Ratio between cross sections with and without Z boson contribution to second 
order in a s . 
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